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ABSTRACT 

Post-Starburst Quasars (PSQs) arc hypothesized to represent a stage in the evolution 
of massive galaxies in which the star formation has been recently quenched due to the 
feedback of the nuclear activity. In this paper our goal is to test this scenario with a 
r~| . resolved stellar population study of the PSQ J0210-0903 ,as well as of its emitting gas 

04 kinematics and excitation. We have used optical Integral Field Spectroscopy obtained 

with the Gemini GMOS instrument at a velocity resolution of « 120kms~^ and spatial 
resolution of « 0.5 kpc. We find that old stars dominate the luminosity (at 4700 A) in 
the inner 0.3 kpc (radius), while beyond this region (at w 0.8 kpc) the stellar population 
is dominated by both intermediate age and young ionizing stars. The gas emission-line 
ratios are typical of Seyfert nuclei in the inner 0.3 kpc, where an outfiow is observed. 
Beyond this region the line ratios are typical of LINERs and may result from the 
combination of diluted radiation from the nucleus and ionization from young stars. The 
gas kinematics show a combination of rotation in the plane of the galaxy and outflows, 
observed with a maximum blueshift of —670 km s^^. We have estimated a mass outflow 
rate in ionized gas in the range 0.3 — 1.1 Mq yr^^ and a kinetic power for the outflow 
of £:out ~ 1-4 - 5.0 X 104°ergs-i « 0.03% - 0.1% x Lboi- This outflow rate is two 
orders of magnitude higher than the nuclear accretion rate of « 8.7 x 10"'^ Mq yr~^, 
thus being the result of mass loading of the nuclear outflow by circumnuclear galactic 
gas. Our observations support an evolutionary scenario in which the feeding of gas to 
the nuclear region has triggered a circumnuclear starburst lOO'sMyr ago, followed by 
the triggering of the nuclear activity, producing the observed gas outflow which may 
have quenched further star formation in the inner 0.3 kpc. 

Key words: galaxies: active - galaxies: kinematics and dynamics - galaxies: stellar 
content - galaxies: starburst 



1 INTRODUCTION 

Our understanding of galaxy evolution has undergone 
a revolution in the past decade. Of particular inter- 
est is the correlation that links the nuclear supermas- 
sive black holes (SMBH) mass to the b ulge stellar ve- 
locity di spersion, the Mbh — cr* re l ation JCebhardt et al.l 
2000allbl: iMerritt fc Ferraresd I2001I: llVeniaine et all |2002|: 



Ferrarese fc Fordll2005l : iGultekin et al.ll2009l : ICraham et all 
2OIII ). In this scenario, Post-Starburst Quasars (PSQs) may 
acquire a great importance, since they seem to represent a 
critical phase in the secular evolution of galaxies that links 
the growth of the stellar bulge and that of the SMBH. 

PSQs are broad-lined AGN that show Balmer jumps 
and high-order Balmer absorption lines from A stars char- 
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acteristic of massive post-starburst populations with age of 
a few hundred Myrs. PSQs are hypothesized to represent 
a stage in the evolution of massive galaxies in which both 
star formation and nuclear activity have been triggered and 
are visible simul t aneou sly before one or the other fades. 
iBrotherton et ahl (| 19991 ) reported the discovery of a spec- 
tacular PSQ that they propose to represent a stage in hier- 
archical galaxy evolution following a merger and enshrouded 
Ultra Luminous Infrared Galaxy (ULIRG) phase. HST im- 
ages have shown that the host galaxy is in fact a post- 
merger remnant, and its spectrum shows evidence for mul- 
tiple stages of star-formation as might be expected through 
such a proc ess jCanalizo et al.ll200a fBrotherton et al.ll2002l : 
ICales et al.ll20m 

There are at least two possibilities to connect the pres- 
ence of a post-starburst population and nuclear activity: (1) 
the flow of gas towards the nucleus first triggers star for- 
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mation in the circumnuclear region and the nuclear activ- 
ity is triggered after hundreds of Myrs; in the meantime, 
the star forinatio n may cease due to exhaustion o f the gas 
IStorchi-Bcrgrnann et alll200ll : lOavies et al.ll2007l 'l: (2) the 
flow of gcis towards the nucleus triggers star formation in 
the circumnuclear region and the nu clear activity, when 
triggered, quenches the star formation ([Granato et al.l 20041: 
Pi Matteo et al.ll2005l : lHopkins et al.ll2006l : ICano-Di'az et alj 
20121 ). In the first scenario, the fueling of the black hole usu- 
ally occurs with a delay of a few 100 Myrs after the begin- 
ning of the star formation - in order to explain the system- 
atic excess of intermediate age stellar population contribu- 
tion in AGN hosts, although there are a few cases in which 
simultaneous on-going star formation and nuclear activity 
are observed (IStorchi-Bergmann et al]|200ll ). In the scenario 
proposed bv lDavies et al. ( 20071 ') the gas fueling the nuclear 
activity originates from winds of AGE stars, which can be 
accreted more efficiently than the ejected mass by O and 
B stars and/or supernova. In these two versions of the first 
scenario, the star formation ceases due to the exhaustion 
of gas, and the post-starburst population is distributed in 
the circumnuclear region extending even to the nucleus. In 
the quenching scenario, the star formation occurring in the 
vicinity of the nucleus is abruptly interrupted due to feed- 
back from the AGN (e.g., nuclear outfiows, radiation and 
jets). Our proposition is that, in this case, the post-starburst 
population should be located in regions suffering the feed- 
back effects - e.g., over regions where a gaseous outflow is 
observed, which can be mapped via the gaseous kinematics. 
With the goal of investigating the nature of the connection 
between the post-starburst stellar population and nuclear 
activity in PSQs, we began a program to map the stellar 
population and the manifestations of nuclear activity in the 
inner few kpc of PSQs using integral field spectroscopy. 

In this paper we present the results obtained from 
observations of the PSQ SDSS J021011. 5-090335. 5 (here- 
after PSQ J02 10-0903), selected for its proximity, at a dis- 
tance of only 170 Mpc (from NEeQ, for Ho = 73.0 ± 
5kmsec~^ Mpc"'^), allowing the study of the spatial dis- 
tribution of its stellar population and gas emission char- 
acteristics. At this distance, the scale at the PSQ is ~ 
0.83 kpc arcsec"^ The PSQ J0210-0903 is also one of the 
brightest PSQs with z <0.1, and its spectrum clearly re- 
veals the presence of the Balmer jump and high-order ab- 
sorption lines of the Balmer series. These absorption fea- 
tures are characteristic of the atmospheres of A stars, which 
dominate the absorption spectra of intermediate age (~ 10* 
yrs) stellar populations (or post-starburst galaxies). The 
PSQ J0210-0903 is hosted by a barre d spiral galaxy with 
Hubble type Sa (|Graham fc Lil \200^ ) and its SDSS Pet- 
rosian absolute magnitude is Mi — —22.32 ±0.50. Although 
this obj ect is not as luminous as typical QSOs (e.g., as 
those in iBrotherton et al.1 12002| : ICales et all 1201 ll ) its ab- 
solute magnitude is cl ose to the lower luminosity limit of 
QSOs (Mi = -22.00; IShen et al.1 120111 ) and its spectrum 
presents similar characteristics to other PSQs. In addition. 



the downsizing issue (jHeckman et al.ll2004l ) makes the PSQs 
with lower redshifts important probes of AGN evolution. 

This paper is organized as follows. In Section [5] we de- 
scribe the observations and reduction processes. In Section 
[3] we present the methodology and results of the stellar pop- 
ulation analysis. In Section |3] we present the methodology 
of analysis of the gas emission and the gas emission-line fiux 
distributions and ratios. In Section [5] we report our results 
for the gas kinematics. In Section[U]we discuss and interpret 
our results as well as discuss the technique of Principal Com- 
ponent Analysis to study the gas kinematics. In sections 16.41 
and 16.51 we present an estimate to the mass of the emitting 
gas and the mass outflow rate, respectively. In Section[7]we 
present a summary of our results as well as our conclusions. 



2 OBSERVATIONS AND DATA REDUCTION 

Two-dimensional spectroscopic data of the PSQ J0210- 
0903 were obtained in 2008 December 7 using the Gem- 
ini Mu lti-Object Spectrograp h Inte gral Field Unit (GMOS- 
IFU) (|Allington-Smith et al.l |2002| ), in the program GN- 
2008B-Q-45. The observations were obtained in two-slit 
mode, covering a field of view (FOV) of 7.0 x 5.2 arcsec^, 
using the B600_G5303 grating and g_G0301 filter, resulting 
in a spectral range of ~ 4400 — 5600 A and a wavelength 
sampling of 0.913 A pixel"^ at a velocity resolution R ~ 
3000 (FWHM~120 kms"^). The total exposure time was 
8x1700 s. The seeing during the observation was approxi- 
mately 0'.'6, corresponding to a spatial resolution at the dis- 
tance of the galaxy of ~ 0.5 kpc. 

Data reduction was accomplished using generic IRAP0 
tasks and specific ones developed for GMOS data in the gem- 
mi. gmos package. The basic reduction steps were trimming, 
bias subtraction, flat-fielding, cosmic rays cleaning, extrac- 
tion of the spectra, sky-subtraction, wavelength and fiux cal- 
ibration, differential atmospheric refraction correction and 
coaddition of different exposures. As no standard stars were 
observed together with the galaxy the fiux calibration was 
only relative. The absolute fiux calibration was obtained 
from the ratio between a spectrum from the Sloan Digital 
Sky Survey (SDSS) database and our spectrum integrated 
within the same aperture (3"). This ratio was then used 
to scale the whole data cube. The final data cube contains 
3640 spectra each corresponding to an angular coverage of 
O'.'l X O'.'l or 83 X 83 pc^ at the distance of the galaxy. Cos- 
mic rays were cleaned from the data before sky-subtraction 
with Laplacian cosm ic ray identification routine LACOSMIC 
(Ivan DokkumllioOll ). The spectra were corrected for redden- 
ing due to the interstellar Galactic medium using the IRAF 
routine noao.onedspec.deredden for the F-band extinction 
Ay = 0.085; its value was calc ulated using the NED extinc- 
tion c alculator, which uses the lSchlegel. Finkbeiner fc Davij 
(|l998h Galactic reddening maps. 

In the top-left panel of Figure [T] we present the Gemini- 
GMOS gr-band acquisition image of the PSQ, which we have 
rotated to the same orientation of the IFU observations. 



1 The NASA/IPAC Extragalactic Database (NED) is operated 
by the Jet Propulsion Laboratory, California Institute of Tech- 
nology, under contract with the National Aeronautics and Space 
Administration. 



^ IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which is operated by the Association of Universities for 
Research in Astronomy, Inc. (AURA) under cooperative agree- 
ment with the National Science Foundation. 
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Figure 1. The top left-hand panel shows a GMOS g-band image of the PSQ J0210-0903 and a box that represents the field-of-vicw 
of the GMOS-IFU observations. The top right-hand panel shows the flux map in H/9 with regions from which the spectra shown bellow 
were identified. Middle panels: spectra at regions 1, 4, 5 and 7 marked at the top right-hand panel showing the emission lines and that 
regions where H/3 has a higher relativity intensity, the continuum is also bluer. Bottom: SDSS spectrum. The scale of the images is 
830 pc arcsec^^ . 



In the top-right panel we present the H/? flux distribution 
within the IFU FOV, including the outline of a set of circu- 
lar regions with O'.'S diameter from which we have integrated 
the spectra used in the stellar population synthesis (see next 
section). In the middle panels of the figure we present rest- 
frame spectra from the nucleus (the AGN, at region 1) to- 
gether with those from two circumnuclear regions (regions 
4 and 5) and from an Hii region, at the more distant region 
7. In the bottom panel the SDSS rest-frame spectrum shows 
the characteristic emission lines of the optical range, as well 
as the absorption lines from intermediate age stars, namely, 
the high-order Balmer lines between 3700 — 4000 A. 



3 STELLAR POPULATION 

We have analyzed the stellar population using both the 
SDSS spectrum — which corresponds to an aperture of 
3" (~ 2.5 kpc at the distance of the galaxy) - covering the 



wavelength range 3770 — 8310 A and our integral field spec- 
troscopy through the use of the technique of spectral syn- 
thesis. 

The Integral Field Spectroscopy could in principle be 
used to obtain information on the spatial distribution of 
the stellar population applying the synthesis to the spectra 
of each pixel. However, the signal-to-noise ratio (SNR) of 
individual spectra was not high enough, thus, in order to 
increase it, we integrated the spectra within the O'.'S aper- 
tures described above and indicated in Figure [1] These re- 
gions were selected as sampling extranuclear regions with 
high enough S/N ratio in the absorption spectra to allow 
performing spectral synthesis and being separated enough 
from the nucleus in order to avoid possible contamination 
from the nuclear light due to the PSF wings. 

In order to pe rform the spectral synthesi s we used th e 
STARLIGHT code (|Cid Fernandes et al.|[200i . |2005| . \200^ ). 
which searches for the linear combination of A'', simple stel- 
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Table 1. Synthesis results. Column 1: identification of spectra. Column 2: the spectral range, where n means narrow range (4630— 
5230 A), m, medium range (4260-5260 A) and ui, wide range (3780-8310 A). Column 3: aperture radius of the spectra. Column 4: 
percent contribution of the AGN featureless continuum to the flux at 4700 A. Column 5-8: population vector in percentage of the 
flux at 4700 A for ages YOUNG, intermediate and OLD. The young population vector was separated in ionizing (xYi'- t ^ 5Myrs) and 
not ionizing {xYnoi- t = 25Myrs) stars. Column 9-11: percent mass fraction for each age. Column 12: stellar mass enclosed in the 
respective aperture in units of solar mass. Column 13: extinction in V-band magnitudes. Column 14: signal/noise ratio (SNR). Column 
15: percent difference between model and observed spectra. 
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lar population (SSP), from a user-defined base, that best 
matches the observed spectrum. Basically, the code fits an 
observed spectrum \ solving the foUowinR equ ation for a 
model spectrum M\ (|Cid Fernandes et al.ll2004 ): 



Mx = Mao | X^a^i&i.ArA 



® G(t). 



(1) 



where Mxo is the synthetic flux at the normalization wave- 
length, X is the population vector, whose components rep- 
resent the fractional contribution of the each SSP to the 
total synthetic flux at Ao, &j,a = Lf'^''{tj,Zj)/Lf^^{tj,Zj) 
is the spectrum of the jth SSP, with age tj and metallicity 
Zj, normalized at \o, rx = 10~'' ''('*^"'*^o' is the redden- 
ing term, and G(w,,cr,) is the Gaussian distribution, cen- 
tred at velocity with dispersion a,, used to model the 
hne-of-sight stellar motions. The reddening term is modeled 
by STARLIGHT due to foreground dust and parametrized 
by the V-band extinction Ay so that all components are 
equally reddened a nd to which we have adopted the G alac- 
tic extinction law of lCardelli. Clayton fc Mathid (| 19891 ) with 
Rv = 3.1. 

The fit of the model to the observe d spectrum is carried 
out m inimizing the following equation (|Cid Fernandes et al.l 
|2004| ): 



M. 



X Wx 



(2) 



where wx is the weight spectrum, defined as the inverse of 
the noise in Ox- Emission lines and spurious features are 
masked out by fixing wx = at the corresponding A. The 
minimum of equation [2] corresponds to the best parameters 
of the model and the search for them is carried out with 
a simulated annealing plus Metropolis scheme. A detailed 
discussion of the Metropolis scheme applied to the stellar 



opul ation synthesis can be found in ICid Fernandes et al.l 
20011 ). 

We constructed the spectral base with the high 
spectral resolution evo lutionary synthesis models of 
iBruzual fc CharlotI (|2002h (BC03), where the SSPs cover 
11 ages, t ^ 1.0 X 10^ 5.0 x 10^ 2.5 x 10^ 1.0 x 10^ 
2.9 x 10*, 6.4x 10^ 9.1 X 10*, 1.4x 10^ 2.5 x lO'', 5 x 10^ and 
1.1 X lO'^" yrs, assuming solar metallicity (Z = 0.02). We 
have us ed the SSP spectra co nstructed from the STELIB 
library (ILc Borgnc et al" 20031 ) , Padova-1994 evolutionary 
tracks and ,Chabrie'ir (| 2003t ) Initial Mass Function (IMF). 
In order to account for the AGN featureless continuum 
(FC) a non-stellar component was also included, represented 
by a power law function (Fi, oc v^^'^). In accordance to 
ICid Fernandes et all (|2004l ) we have binned the contribution 
of the SSPs into three age ranges: YOUNG {xy = 2;Yi -f ai^Ynoi, 
where xyi (tS55 Myrs) represents the ionizing population and 
aiYnoi {t = 25 Myrs) represents the non- ionizing one), inter- 
mediate {xi to 100 Myr ^ t ^ 2 . 5 Gyr ) and OLD {xq to 
t > 2.5 Gvr). [Cid Fernandes et all (|2005l ) have shown that 
the star formation history of a galaxy may be very well 
recovered by this condensed population vector, as the in- 
dividual contributions of each SSP are very uncertain. We 
have used this conservative spectral base, because we do not 
have sufficient constraints in our spectra to allow the use a 
broader spectral base. This proc edure has been extensively 
used in previous studies (e.g., ICid Fernandes et al.l |20M . 
l2005h and the investigation of the impact of the correspond- 
ing assumptions in the fit of the observed spectra is beyond 
the scope of the present paper. Our goal here is to investi- 
gate the relative contribution of broad age components as a 
function of location in the observed field. 

The results of the synthesis are shown in Table [1] and 
Figure [21 they are presented in two sets, one for the SDSS 
spectrum and another for the IFU integrated spectra, which 
are discussed in the next two sections. 
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Figure 2. In this figure we present tlie synthesis results for the SDSS spectrum, for the nucleus (IFU 1), for the circumnuclear regions 
(IFU 4, 5 and 6) and for the Hll region located slightly farther from the nucleus (IFU 7). The black dashed-dotted line at the top of 
each panel shows the observed spectrum. The red solid line shows the fitted model, while the dotted line shows the masked regions of 
the spectrum (emission lines or spurious features). At the bottom we show the spectra corresponding to each age bin (young - blue; 
INTERMEDIATE AGE - green; OLD - red) scaled to its percent contribution to the total continuum light at 4700 A. 



3.1 Synthesis of the SDSS spectrum 



The SDSS spectrum shows clear signatures of intermediate 
age stars which contribute approximately 13% of the total 
luminosity at 4700 A according to the synthesis (first line of 



Table[T]). The old population contributes with 36% while the 
young and FC contribute with 18% and 32% respectively. 

Our IFU data cover a limited spectral range, from 4260 
to 5230 A, thus we investigate the effect on the synthesis 
of this restricted range using the SDSS spectrum. We have 
labeled as m the spectrum covering this range in Table 1. 
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Figure 3. Integrated fluxes in the [O III ] and H/3 emission line in logarithmic units of erg cm~^ , obtained from GH fitting using the 
Levenberg-Marquardt least squares. The errors in the [Olll] flux map are ft! 4% within the inner arcsecond, reaching up to 13% at the 
borders. In the H/3 flux map the errors within the inner arcsecond are ~ 11%, increasing to Si25% at the borders of the mapped field. 



For this spectral range the synthesis results show an increase 
in the contribution of the old population, which is 41%, as 
well as of the young and FC components, which contribute 
with 33% and 27%, respectively, while the contribution of 
intermediate age stars disappears. Outside the nucleus, the 
blue end of the spectra are quite noisy, and we had to re- 
strict the spectral range to the interval 4630-5230 A, which 
we have labeled n in the Table [T] In order to verify the effect 
of restricting the spectral range even further, we performed 
the synthesis of the SDSS spectrum also for this range. The 
result is a further increase of the old stellar population con- 
tribution to 59%, and a decrease of the young population 
and the FC component, while the intermediate age popula- 
tion continues to be absent. 

In summary, using the wider spectral range we find that 
the intermediate age population contributes with approxi- 
mately 13% of the flux at 4700 A within the inner 3 arcsec 
of our PSQ. However, when we restrict the spectral range 
to the intervals labeled as m and n we lose the sensitivity to 
the intermediate age stellar population, most probably due 
to the fact that the main signatures of intermediate age stars 
(high order Balmer absorption lines) have been removed. 

3.2 Synthesis of the IFU spectra 

We have first constructed an "integrated" spectrum by 
adding the contribution of all spectra within the SDSS aper- 
ture (3"). We have then performed the spectral synthesis in 
order to check the consistency of our results. These results 
are shown in the fourth line of Table [T] for the m wavelength 
interval, and they approximately agree with those for the 
SDSS, but presenting a diflerence of ~10% in the population 
vectors xo and xfc- According to the adev parameter, the 
percent difference between the model and the observed spec- 
trum is ~2% for both the SDSS and the IFU spectra, which 
corresponds to a good fit. However, because the stellar pop- 
ulation synthesis method does not provide error estimates 
for the population vector, we can use the difference of «10% 
for the contribution of each component to the total hght at 
4700 A as an error estimate for the population vector. In 
the 5**" to the 12*'' line of Table[T]we show the results of the 



synthesis for the set of circular regions described in section 
[2] corresponding to apertures of O'.'S diameter. As discussed 
above, our IFU data spectra cover a limited spectral range, 
namely the ranges labeled as m and n. Only at the nucleus 
(region 1) we had enough SNR to cover the m range, while 
in all other regions (2 to 7) we use the range n. 

In the nucleus (region 1), we have a predominance of 
the old population, contributing with 59%, while the inter- 
mediate age population contributes with 27% and the young 
-I- FC contribute with 14%. For the circumnuclear regions 
there is an increase in the contribution of the young popu- 
lations, as we can see in the results for regions 2, 3, 4 and 
5. The contribution of the old population decreases from 
59% — at the nucleus — to 54%, 48%, 29% and 37% re- 
spectively, showing that the old population has the largest 
contribution at the nucleus. At region 6, only the intermedi- 
ate age component is present. The young and intermediate 
age stellar populations are thus distributed preferentially in 
the circumnuclear region, at about 0.8 kpc from the central 
AGN. 

Figure [2] shows the results of the synthesis for the SDSS 
spectrum and for five of the IFU integrated spectra of the 
circumnuclear region of the galaxy. The black dashed-dotted 
line, at the top of each panel, shows the observed spectrum. 
The red solid line shows the fitted model while the dotted 
line shows the masked regions of the spectrum. We have 
masked out emission lines and spurious features, such as the 
CCD gap. In each panel, at the bottom, we show the spectra 
corresponding to each age bin (identified as YOUNG, inter- 
mediate AGE and old) scaled to its percent contribution to 
the total continuum light at 4700 A. 



4 EMITTING GAS 

4.1 Measurement of the emission lines 

The kinematics and the flux distribution of the emitting gas 
in the central region of active galaxies have typically been 
measured using a single Gaussian function for each emission 
line proflle. Nevertheless, a close inspection of the emission 
lines from the Narrow Line Region (NLR) of AGNs (e.g.. 
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iRiffel et al] l2009al . l2010al ') shows that the profiles present 
deviations from a purely Gaussian profile, having "blue" or 
"red" wings, that are usually modeled by fitting a set of 
two or more Gaussians. In order to take into account the 
presence of these asymmetries, we used instead of Gaussians 
the Gauss-Hermite series (hereafter, GH), with which we 
have fitted the emission-line profiles of H/3 and [O iii jfj- 

The GH series allows us to fit the wings of the emis- 
sion lines via the /is and /i4 moments, in addition to the 
centroid velocities and the v elocit y dispersion. Following 
Ivan der Marel fc Franxl (|l993l 'l and ICappellari fc Emselleml 



[Olll]/H^ 



( 20041 ) ■ the GH series can be written as 



A- Ao 



(3) 



where 7 is the amplitude of the GH series, Ao is the central 
wavelength, a is the velocity dispersion, hj are the GH mo- 
ments and Hj are the Hermite polynomials and the function 
a(w) is the standard Gaussian, given by 

1 _i£ 



(4) 



For any choice of free parameters of the Gaussian (7, Ao 
and o") there is a GH series described by the equation ((3]) 
that fits the emission line profile. A series that converge 
rapidly to the desired solution and that is very similar to the 
standard Gaussian function can be obtained by summing up 
to N=4 and t aking the assumption that feo = Ho{w) = 1 and 
/ii = /i2 = Jvan der Marel fc Franxll 19931 ') . In this case, the 
GH series given by equation ((3| can be approximated by 

a{w) 



GH{w)=j- 



[1 + hsHaiw) + hiH^iw)] 



(5) 



where H3{w) and H4,(w) are the 3rd and 4th order Hermite 
polynomials. 

The /13 moment measures asymmetric deviations of a 
Gaussian profile, such as "blue" (h^ < 0) or "red" (/13 > 0) 
wings. The hi measures the symmetric deviations, where a 
wider and flatter profile than a Gaussian has /i4 < and 
narrower and more peaked profile than a Gaussian has hi > 
0. In general, the use of the GH series allows a better fit to 
any emission-line profile. 

The H/3 and [Oiii] emission- line profiles were measured 
using an IDL routine that solves for the best solution of pa- 
rameters of the GH series using the nonlinear least-squares 
Levenberg-Marquardt method . The ID L routine used is sim- 
ilar to that detailed in iRiffell (l2010bh which makes use of 
the routine MPFIlH (|Markwardtll2009l ). which performs the 
minimization. We have performed the fit to the spectrum of 
each pixel of the data cube in order to obtain the spatial 
distribution of the emission-line fluxes, velocity dispersions, 
centroid velocities and GH moments /13 and /14. 

4.2 Error estimates 

Error estimates were obtained using Monte Carlo simula- 
tions adding to the original spectrum an artificial Gaussian 



^ We denote the [Olll]A5007 emission line simply by [Olll] 
* MPFIT provides a robust and relatively rapid way 
to run the least squares method and is available in 
|http: //purl.com/net/mpfit| as part of Markwardt IDL Library 
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Figure 4. Distribution of the emission-line ratio [O III ]/H/3. The 
regions from which spectra have been extracted for the study of 
the stellar population are identified by the circles. 



noise normally-distributed, whose mean is equal to zero and 
the standard deviation equal to one. We have performed the 
Monte Carlo simulations with 200 realizations computing 
the best-fitting parameters of each realization so that we 
have the standard deviation of the best model in the end 
of the simulation. In principle, the larger the number of re- 
alizations, more reliable are the results, but in our case we 
have verified that with approximately 200 realizations the 
standard deviations are already of the same order as those 
obtained with 500 to 1000 realizations. In this way, we have 
obtained the error estimates for the free parameters of the 
GH function: the centroid velocity, velocity dispersion, the 
moments /13 and /14 and the flux of the emission line. The 
error values in each measured parameter are given together 
with the description of the corresponding maps in the fol- 
lowing sections. 

4.3 Emission line flux distributions 

In Figure [3] we present the flux distributions in the H/3 
and [Oiii]A5007 emission lines, both given in logarithmic 
scale and erg cm^^ s~^ units. The location of the maximum 
brightness of the [O iii ] flux distribution coincides — consid- 
ering the spatial resolution of our data — with the position 
of the peak of the continuum. This position was determined 
as the centroid of the flux distribution in the continuum, ob- 
tained by collapsing the data cube between two continuum 
wavelengths. We have adopted this position as the galaxy 
nucleus and its location in the figures is identified either by 
two perpendicular dotted lines or a cross. 

The H/3 map shows a fiux distribution elongated to- 
wards the South of the nucleus with a detached patch to the 
North. In the lower right corner there is another patch of 
emission which we attribute to an Hii region located within 
what appears to be a spiral arm of the galaxy. The [O iii ] 
flux distribution is also elongated along the North-South, 
but along the East-West is broader and m ore symmetri- 
cally d istributed around the nucleus than H/3. lGraham fc Lil 
(|2009D . using SDSS images, have reported the presence of a 
bar along the North-South direction, where we have reported 
the elongations in [Oiii] and H/3 flux distributions. The er- 
rors in the [O iii ] emission line fluxes are ~ 4% within the 
inner arcsecond, increasing to « 13% at the borders of the 
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Figure 5. Left: Centroid velocity maps of [Olll] and H/3 in units of kms^^. Right: Velocity dispersion maps. Both have been obtained 
from the fit of GH series to the emission-line profiles. The mean uncertainties are between ^6kms^^ and 15kms^^ for both the velocity 
fields and velocity dispersion maps for [Olll]. For H/3 the mean uncertainties are between RilOkms^^ and 25kms^^. 



mapped field. In the case of H/3 the errors within the inner 
one arcsecond are « 11% while at the borders of the mapped 
field they reach « 25%. 

In order to map the gas excitation, we show in Figure [4] 
the ratio between the [Oiii] and II/3 fiux maps, whose dis- 
tribution presents the highest ratios — and thus the highest 
excitation — at the nucleus and up to « 1" to the North, 
roughly within the circle (with diameter of O'.'S) labeled by 
the number 1. In this region the [Oiii]/H/3 ratio reaches a 
value of ^ 20, which decreases steeply outwards, down to 
values as low as ~ 3, which is typical of LINERs and may 
include also contribution of ionization by young stars. The 
lowest values of [Oiii]/H/? are found at the regions labeled 
4, 5 and 6, in regions where the stellar population synthe- 
sis results showed the largest contribution of young stellar 
populations. 



5 GAS KINEMATICS 

As described in the previous section, we have obtained the 
gas kinematics — centroid velocities, velocity dispersions 



and GH moments and /14 — using the GH series to fit 
the emission lines H/3 and [Oiii]. In Figure [S] we present, 
to the left, the centroid velocity maps and, to the right, the 
corresponding velocity dispersion maps. The velocity maps 
are masked according to the error estimates, which are on 
average in the range 10 — 20kms^^. Pixels for which the 
error values are larger than these were masked. Centroid ve- 
locities are shown relative to the systemic velocity of the 
galaxy, 12480±20 kms~^, whose value was measured from 
the velocity channel maps in the H/? emission-line, as ex- 
plained in section 15.11 

The [O III ] velocity field (top left panel of Fig. (5)1 shows 
positive values towards the South and negative values to- 
wards the North, suggesting a rotation pattern with the line 
of nodes approximately along the North-South direction and 
a velocity amplitude of « 150 km s~ ^ . Alternatively, this pat- 
tern could also be due to a bi-polar outflow, frequently ob- 
served in the Narrow-Line Region of AGN. Considering that 
the velocity pattern is distorted relative to that correspond- 
ing to simple rotation and it's not centered at the nucleus, it 
is more likely that both rotation and outflows are present, as 
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Figure 6. Distribution of the GH moments (left panels) and /14 (right panels) of [Olll] (top panels) and H/3 (bottom panels) emission 
lines. While the moment measures asymmetric deviations of a Gaussian profile, such as "blue" (/13 < 0) or "red" {h-s > 0) wings, the 
/14 measures the symmetric deviations, where a wider and flatter profile than a Gaussian has /14 < and a narrower and more peaked 
profile than a Gaussian has /i4 > 0. The average errors for the /13 and /i4 in [Olll] are ^0.02 within the inner arcsecond and 0.08 at 
the borders. The errors for /13 and hi in the H/3 arc larger than for [Olll], increasing from 0.04 within the inner arcsecond to 0.2 at the 
borders. 



evidenced also by the blueshifts around the nucleus. The av- 
erage errors in the velocity values increase from ~6kms^^ 
within the inner arcsecond to « 15kms~^ at the borders of 
the mapped field. 

The [O III ] velocity dispersion map (top right panel of 
Fig. [5} presents a range of values from »il30kms~^ to 
RiSOOkms"'^. Within a radius of 0'.'5 from the nucleus we 
observe the lowest cr's, of «130kms~^. This region is sur- 
rounded to the West by a semi-circle where the a values 
are slightly higher, « 175 km s"'^, while to the SE, East and 
NE there is another semi-circle where the values reach even 
higher values of « 250-300 kms~^. The errors and their dis- 
tribution are similar to the ones for the [O iii ] velocity field. 

The H/3 velocity field (bottom left panel of Fig. [5]) also 
presents a pattern that suggests rotation plus outflows with 
similar orientation to that observed in [O iii ] , but reaching a 
higher velocity amplitude of ~200kms~^. As in the case of 
[Oiii], we also observe a rotational pattern with kinematic 
center displaced from the continuum peak, showing blueshift 
at the nucleus. The average errors in the H/3 velocity field 
range from lOkms"^ for the inner arcsecond to 25kms~^ 
at the borders of the mapped field. 

The H/3 velocity dispersion map (bottom right panel of 



Fig. [5]) shows a smaller region of 0'.'2 near the the nucleus 
where the a has values of ^ISOkms^'^ surrounded by a 
semi-circle to the West where the a is larger than in the 
case of [Oiii] (reaching ~200kms~^). Outwards, the values 
decrease to ~100kms~^. Although the H/3 emission could 
not be measured as far as that of [Oiii] towards the East, 
the small region where the cr(H/3) could be measured also 
shows an increase there, similar to what we have observed 
for o-([OIII]). 

In the top panels of Fig. [6] we show the /13 (left) and 
hi (right) GH moments of the [Oiii] emission-line profiles. 
The /i3 moment, over most of the IFU field, presents values 
between and —0.1, indicating that there isn't significant 
asymmetric deviations with relation to Gaussian profiles. 
However, it reaches values of up to —0.3 to the NE, indicat- 
ing the presence of blue wings in this region, giving addi- 
tional support to the presence of outflows. The /14 moment 
shows only positive values, meaning that the profiles are 
narrower and more "pointy" than Gaussians. In the central 
region the moment has values around 0.2, increasing to 
0.3 to the NE and decreasing to about 0.1 to the SW. The 
average errors obtained from the Monte Carlo simulations 
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Figure 7. Channel maps obtained from slices - of width with 55 kms~^ - along the emission-Une profiles of H/3. The numbers in the 
top left of each panel are the central velocities of each slice in kms~^. Fluxes are shown in logarithmic scale and in units of ergcm"^ s~^. 
Spatial positions are in arcseconds. The solid black line represents the major axis of the galaxy. 



for the hz and /14 are «0.02 within the inner arcsecond, 
increasing to 0.08 at the borders of the mapped field. 

In the bottom panels of Fig. [6] we show the /13 (left) 
and /i4 (right) GH moments of the H/3 emission- line profiles. 
The hz moments present values of about —0.3 around the 
nucleus, suggesting the presence of outflows towards us in 
this region, whereas in other regions the values are close to 
zero. The /i4 moment presents an irregular distribution, with 
values not deviating much from zero, indicating there are no 
significant symmetric deviations from a Gaussian profile for 
the H/3 emission line. The average errors in hz and /i4 are 
larger than for [Oiii], increasing from 0.04 within the inner 
arcsec to 0.2 towards the borders of the mapped field. 



5.1 Channel maps 

We have obtained channel maps by integrating the flux dis- 
tribution in different slices of velocity along the [O iii ] and 
emission-line profiles. We have integrated the flux in each 
velocity channel after subtraction of the continuum contri- 
bution from both sides of the profile. This technique allows 
us to map the velocity field of the emitting gas throughout 
the emission line and not only at the central wavelength (like 
the velocity maps shown in Fig. [5]). 

In Figure [T] we present the channel maps of the H/3 



emission line for f6 velocities. The systemic velocity has 
been subtracted from the maps and its value was obtained as 
described below. The highest blueshifts, of « — SSOkms"'^, 
are observed at the nucleus. As the velocities become less 
negative (channel —273 km s^^), it is possible to observe 
emission outside the nucleus towards the North. For veloc- 
ities around zero, we can observe emission distributed over 
most parts of the field, including in the spiral arm to the 
West of the nucleus and in a similar structure to the East. 
As the velocities become more positive, the emission distri- 
bution moves to the South side. The highest positive ve- 
locities reach 290kms~^ to the South of the nucleus. The 
kinematics observed in the velocity channels are consistent 
with a rotation pattern with line of nodes approximately 
along North-South and an outflow observed in negative ve- 
locities within ~f arcsec from the nucleus. As the nuclear 
outflow disturbs the rotational pattern in the nuclear re- 
gion in the centroid velocity maps (Fig. [S]), we have used 
the H/3 velocity channels to estimate the systemic velocity, 
observing that the flux distribution pattern becomes sym- 
metric relative to the nucleus along the galaxy major axis 
for a systemic velocity of 12480 ± 20kms~'^. A comparison 
between the panels at —273 and -1-290 km s^^, for example, 
(approximately symmetric relative to the zero velocity) in- 
deed shows that the extranuclear emission is approximately 
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Figure 8. Channel maps obtained from slices - of width 55 km s~^ — along the emission-line profiles of [O III ] . The numbers in the top 
left of each panel are the central velocities of each slice in kms~^. Fluxes are shown in logarithmic scale and in units of ergcm~^s~^. 
Spatial positions are in arcseconds. The solid black line represents the major axis of the galaxy. 



symmetric relative to the nucleus (within ~20kms~^). In 
addition, most of the extranuclear emission along the mi- 
nor axis is found at zero velocity, supporting our adopted 
systemic velocity. 



In Figure [51 we present the channel maps of the [O iii ] 
emission line in 25 different velocity channels, with an in- 
crease in velocity of « 55 kms~^ from one panel to another. 
The highest negative velocities, that reach — 763kms~^ — 
therefore higher than the ones observed for H/3 — are also 
observed at the nucleus, but extending circa I'.'O towards 
Northeast (NE) of the nucleus. As the velocities become 
less negative, the emission becomes more extended toward 
the North. For velocities around zero we can observe the 
emission extending from North to South. As the velocities 
become more positive, the emission moves towards South, 
in the same way as observed for H/3 (up to 275 kms~^). For 
[O III ] even higher positive velocities are observed, however 
the peak values remain centered on the nucleus. 



6 DISCUSSION 



6.1 Gas Excitation 



We have investigated the gaseous excitation by mapping the 
ratio of [Oiii]/H/3 and by measuring its value in the in- 
tegrated spectra of the circumnuclear region, as shown in 
column 6 of Table [5] Additionally, we have separated the 
synthesis results for the young stellar population in terms of 
ionizing (xyi) and not ionizing population (a;Ynoi), given in 
percent contribution in columns 5-6 of Table [T] 

Using as reference the B PT d ia gnostic diagrams 
dBaldwin. Phillips fc TerlevichI Il98ll: Kauflmann et all 



l2003l : iKewlev et al.ll2006l : IStasiriska et al.ll2006l ') we conclude 
that the emission-line ratios within the inner 0.3 kpc 
(radius) are typical of Seyfert nuclei with a mean value 
of [Oiii]/H/3 =17.9, but reaching up to ~25 at locations 
closest to the nucleus. In the circumnuclear region smaller 
values are observed, decreasing from 8.1 (at region 2) to 0.2 
in the detached region 7. At these locations (regions 2-7), 
the synthesis shows that the young population is essentially 
comprised by ages younger than 5 Myrs. The BPT diagrams 
reveals additionally LINER-like emission-line ratios, which 
we have interpreted as due to excitation by the combined 
contribution of radiation from the AGN and hot stars. 
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Figure 9. Tomograms 1, 2 and 4 of the PCA analysis (top) and their respective eigenspectra (bottom). The eigenspectra show the 
correlations and anti-correlations in the data while the tomograms map them spatially. Positive and negative features in the eigenspectra 
are associated with regions in the tomograms whose color scales show respectively positive and negative features. The wavelength scale 
in the eigenspectra are given in the galaxy's reference frame. 



6.2 Gas Kinematics 

The gaseous velocity fields of [Oiii] and H/3 presented in 
Fig. [5]botli show similar rotation patterns, but with the ro- 
tation displaced from the continuum peak position by circa 
I'.'O to the South. Such velocity field suggests that a possi- 
ble rotation pattern around the nucleus is disturbed by a 
kinematic component not due to the gravitational potential 
of the bulge. In this region the emitting gas is blueshifted, 
with the velocity centroid reaching values of the order of 
— 200kms~^. The velocities observed in the wings of the 
emission-line profiles are, nevertheless, much higher, reach- 
ing — 760kms~^ for [Oiii] and — 670kms~^ for H/3, as ob- 
served in the channel maps of Figs. [8] and [T] These blueshifts 
support the presence of a nuclear outflow, in agreement also 
with the high velocity dispersions observed in this region 
(right panels of Fig. [5|. Considering both the centroid ve- 
locities and the channel maps, the kinematic of the emitting 
gas can be understood as the combination of gas in rota- 
tion in the plane of the galaxy combined with an outflow 
component. 

The orientation of the rotation axis is approximately 
along the North-South direction, although the true orienta- 
tion is not clear due the presence of the outflow in the central 
region. The rotation in the H/3 emission is seen in the chan- 
nel maps between —273 and -|-290kms~^, with blueshifts 
observed along the major axis to the North of the nucleus. 
These blueshifts decrease gradually as the emission shifts 



towards the South side of the major axis and the velocity 
increases to redshifts. 

Additional support for the rotation in the [O iii ] cen- 
troid velocity field is given by the Principal Component 
Analysis of the data cube, which is described in the next 
section. The Principal Component 4 (PC4) together with 
its tomogram seems to be consistent with the presence of a 
rotational component, as shown in the channel maps of the 
[O III ] emission-line. 

6.3 PCA tomography 

Astronomical observations with IFU's provide a large 
amount of information, with two spatial dimensions and one 
spectral dimension, so that there are tens of millions of pix- 
els to be analyzed. An efficient method to analyze all this 
information in an optimized way is the Princip al Component 
Analy sis (PCA), which has been adapted bv ISteiner et al.l 
(|2009h for application in IFU data. The PCA technique has 
been previously employe d in astronomy fo r morphological 
classification of galaxies (|Lahav et al.lll996l ) and for distin - 
guishing radio-loud from radio-quiet QSOs (|Borosonll2003 '). 

In order to extract additional information from our data 
and its complex gas kinem atics, we have employed the PCA 
technique as described bv lSteiner et al.l (|2009l ). to the cali- 
brated data cube. The PCA allows us to separate the infor- 
mation initially contained in a system of correlated coordi- 
nates into a new system of uncorrelated coordinates, whose 
components are ordered in accordance with decreasing val- 
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ues of the variance of the principal components. The new 
coordinates are called "eigenspectra" — i.e., eigenvectors 
given as a function of wavelength — so that their projec- 
tion on the data produces images which are called "tomo- 
grams". While the eigenspectra show the correlations and 
anti-correlations in the data, the tomograms will map them 
spatially. The tomograms are like slices of data in the space 
of eigenspectra, whose interpretation can reveal correlations 
or anti-correlations between spectral and spatial character- 
istics of the data. 

In the panels of Figure [9] we present the results of the 
PCA applied to our data cube, where we show the princi- 
pal components PCI, PC2 and PC4 (which show the largest 
variance of the data), represented through the correspond- 
ing tomograms and eigenspectra. We are not presenting the 
PCS, which has the third largest variance, because it con- 
tains mainly noise, mainly due to the difference in response 
between the blue and red slit detectors of the IFU. The other 
components do not present a significant variance. Correla- 
tions in the emission lines appear as all of them having pos- 
itive values in the eigenspectrum (e.g., bottom left panel of 
Fig. [9}. Anti-correlations appear as one line having positive 
value and the anti-correlated one having negative value (e.g., 
H/3 negative and [O iii ] positive in the bottom middle panel 
of Fig.[9l). 

In the left panels of Figure [9] we show tomogram 1 
and its eigenspectrum (PCI), which represents the princi- 
pal component of the PCA with circa 93.4% of the total 
variance. It shows a positive correlation between the emis- 
sion lines, indicating that they are all generated in the same 
spatial region — the brighter region of the tomogram. Since 
the emission lines that are present in the eigenspectrum have 
[O iii]/II/? ratio typical of ionization by AGN, as well as con- 
tribution of the stellar population, we conclude that PCI re- 
veals the dominance of the AGN and the stellar population 
over the entire central region. 

The central panels of Figure [9] show eigenspectrum and 
tomogram 2, for which there exists an inverse correlation 
between [O iii ] and both the H/3 emission and the stellar 
component, whose features are inverted and below zero in 
the eigenspectrum. The tomogram indeed shows that the 
corresponding emission originates in different regions: [O iii] 
dominates in the inner « O'.'S (positive values), while the H/3 
and stellar population dominate in the surrounding region 
(negative values). The interpretation that we propose for 
PC2 is that although PCI shows that most of the H/3 emis- 
sion is due to the AGN, there is some H/3 emission which 
has a different origin. This emission is thus not generated 
by the AGN, but possibly by young stars that ionize the 
gas, indicating the presence of recent star formation around 
the nucleus. This conclusion is in agreement with the results 
of the stellar population synthesis, which show, in fact, the 
presence of young stars in these regions. 

In the right panel of Figure we show PC4, which rep- 
resents 0.09% of the variance. The eigenspectrum 4 shows 
that the [O iii ] emission has part of the profile in redshift and 
part in blueshift, and that they are anti-correlated, since the 
blueshift profile has negative values and the redshift profile 
has positive values; i.e. they originate in different regions, 
as seen in tomogram 4. The part in redshift originates in 
a compact structure to the Southwest of the nucleus which 
curves and extends to the South, while the blueshifted part 



originates in a compact region to the Northeast of the nu- 
cleus. This structure may be identified with the rotation 
component seen in the [Oiii] channel maps. 

6.4 Mass of the emitting gas 

We have estimated the total mass of ionized gas of the inner 
region (excluding the detached spiral arm) as: 

M^Ven^rrip, (6) 

where rio is the electron density, mp is the proton mass, e is 
the filling factor and V is the volume of the emitting region. 

The filling factor can be calculated using the H/3 lumi- 
nosity, as follows. The emiss ivity of the H/3 line is given by 
(lOsterbrock fc Ferlandl [200^ 1: 



off hVYip 

■ ncnpOHi^^— 



(7) 



where rip is the proton density, is the H/3 effective re- 

combination coefficient and hi/n/s is the corresponding en- 
ergy. The total luminosity in H/3 is given by integrating the 
emissivity over the entire emitting volume and over all di- 
rections, assuming a pure hydrogen gas which is completely 
ionized (no = rip): 



L(H/3) 



jufi dfldV = f .24 X 10" 



erg s 



(8) 

where we have adopted the a'^^J for a temperature of lO** K. 
From the above equation, we obtain the product Ve: 

-,59 -^-41 (H/J) -3 



Ve = 8.1 X 10" 



(9) 



where L41 (H/3) is the H/3 luminosity, in units of lO'*^ ergs s~^ 
and ria is the electron density in units of 10'' cm""^. Thus, 
we can estimate the mass of the emitting region by inserting 
the expression [9] for Ve above in|6] obtaining: 

,5^41 (H/3) 



/W" « 7 X 10" 



Mr, 



(10) 



given in units of solar masses (|Petersonlll997l ). 

The H/3 luminosity was calculated from the inte- 
grated H/3 flux F(H/J) of the inner region corrected 
for the reddening obtained from the SDSS spectrum 
(C(H/3) = 1.38 ± 0.05), using the reddening law of 
ICardelli. Clavton fc MathisI (|l989t ) adopting the theoretical 
ratio of 3.0, corresponding to case B recombination 
Osterbrock fc Ferlandl 120061 ) . For the assumed distance of 
d =170Mpc, we obtain L(H^) = 47rd^f (H/3)10'^'"''' = 
1.79 ±0.21 X lO^^ergss-^ Using the [S 11] A6716/A6731 ra- 
tio, obtained from the SDSS spectrum, the electron den- 
sity ns was obtained by solving numerically the equilibrium 
equation for a 5-level atom using the IRAF routine STS- 



DAS . ANALYSIS . NEBULA R . TEMDEN (jPe Robertis et al] Il987l : 
IShaw fc DufouJll994ll The assumed electron temperature 
was 16000 K (|Petersonlll997^ ■ resulting in na = 0.150±0.02. 
With these assumptions, we obtain a mass of ionized gas of 
8.34 ± 0.92 X W^Mq. 



6.5 Mass outflow rate 

We have concluded from the gas kinematics that the 
blueshifted emission in the nuclear region is due to out- 
flowing gas. Although the geometry of the outflow is 
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not c lear from our maps , prev i ous kinematic studies 



fe.g.. ICrenshaw fc Kraemed I2OO0I : iDas et alj |2005| . I2OO6I : 



IStorchi-Bergmann et al.ll2010D of other active galaxies show 
that AGN outflows present a conical geometry. On the ba- 
sis of these previous studies, we assume that the outflow in 
PSQ J02f0-0903 also has a conical geometry, with the cone 
axis directed approximately towards us. We thus consider 
that the outflowing gas is crossing the base of a cone whose 
radius we estimate from the extent of the blueshifted region 
around the nucleus seen in the channel maps. We have mea- 
sured the ionized gas mass outflow rate at negative velocities 
through a circular cross section with radius r = 0'.'6 around 
the nucleus, assuming that the height of the cone is equal 
to the diamet er of its base. The mass outflow rate can be 
calculated as (jRiffel fc Storchi- Bergmann|[2011al lbh: 



Mout — mpUeVo^itAe , 



(11) 



where A = nr « 7.42 x 10 cm is the area of the circu- 
lar cross section, nip = 1.67 x 10~^*g is the proton mass, 
rie is the electron density and e is the filling factor. The 
filling factor can be obtained from equation |9] above tak- 
ing into account only the luminosity of the gas in outfiow. 
We have considered that the outflowing gas corresponds to 
that observed in the velocity channels between —667 and 
-104kms"\ The total flux was obtained by adding the flux 
of the channels centered between —667 and — 329kms~^ 
and by adding 50% of the flux of the channels between 
—273 and — 104kms^^, for which the remaining emission 
- the other 50% - is assumed, as a best guess, to originate 
in the plane of the galaxy. The total fiux in these chan- 
nels is F(H/3) = 1.68 X 10"^* ergcm'^ s"^ and L(H^) = 



47rd^F(H/3) = 5.81 x 10' 



ergcm"'^ 

*° ergs~^. From the adopted geom- 
etry, the volume of the cone whose height h and radius r are 
such that /i = 2 r = 996 pc is = 7.6 X lO*'^ cm^. The veloc- 
ity adopted for the outfiow is the average of the velocities 
of each channel above weighted by the corresponding fiux, 
which results Vout = — 340kms~^. Adopting Ue = 150 cm""^ 
(obtained from the [Sii] line ratio from the SDSS spectrum), 
e = 0.0027, we obtain Mout = 1.1 yr~^. For a somewhat 
larger gas density, of Ue = 500 cm~^ (as the aperture of the 
SDSS is too large, and the inner gas density being probably 
larg er), we obtain e = 0.00025 a n d Mout = 0.3 Mq yr"^ 

iRiffel fc Storchi- BergmannI (|2011bh have obtained 
Mout = 6 M0 yr~^ for Mrk 1157. They have a lso revised pre- 
vious results for another 6 active galaxies in iBarbosa et al.l 
l|2009l) and in lRiffel et al] (|2009bh . obtaining values for the 
mass outfiow rate in range of Mout ~ 0.1 — 6M0yr~^. 
Similar values, of ~ 2.0M(7)y r~^, have been reported by 
Storchi-Bergniann et al] (|2010l ') for NGC 4151 as well as by 
Veifleux et al.l (120051 ) . with Mout ~ 0.1 - lOMoyr-^ for 
a set of luminous active galaxies. The mass outflow rate 
that we have estimated for PSQ J0210-0903 is thus in good 
agreement with previous estimates for other AGNs. 

Now we can compare the outflow mass rate with the 
accretion r ate necessary to feed the AGN, which can be cal- 
culated as (jPetersonll 19971 ): 



Lhoi 
rj 



1.8 X 10" 



— ) Moyr- 
77 



(12) 



where Lboi is the nuclear bolometric luminosity, c is the light 
speed and r] is the efficiency of conversion of the rest mass 
energy of the accreted material into radiation power and L44 



Table 2. Mean stellar ages. Column 1: identification of spec- 
tra. Column 2: the spectral range, as in Table [T] Column 3: 
aperture radius in arcsec. Column 4 and 5: mean stellar age 
weighted by the light and mass fraction, respectively. Column 
6: Emission line ratio [Olll]/H/3 for corresponding spectrum. 



Pes. 


AA 


ap. 


<logt, >L 


<logt, >M 


[Oni]/H/3 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 






(arcsec) 


(yr) 


(yr) 




SDSS 




1.5 


8.46 


9.94 


8.46 


IFU 1 


m 


1.5 


8.73 


9.84 


11.15 


IFU 1 


m 


0.4 


9.48 


9.96 


17.9 


IFU 2 


n 


0.4 


8.23 


10.03 


8.1 


IFU 3 


n 


0.4 


7.94 


10.03 


7.2 


IFU 4 


n 


0.4 


7.48 


9.99 


4.6 


IFU 5 


n 


0.4 


8.22 


9.98 


3.4 


IFU 6 


n 


0.4 


8.93 


8.95 


2.6 


IFU 7 


n 


0.4 


6.00 


6.00 


0.2 



is bolometric luminosity in units o f 10**ergs~^. Following 
iRiffel fc Storchi-Bergmamil (|2011al lbh. I/boi can be approxi- 
mated as « lOOL(Ha), where L(Hq) is the Ha nuclear lu- 
minosity. Using the interstellar extinction coefficient C(H/3), 
the intrinsic ratio ffSsl ^ 3.0 (jPsterbrock fc Ferland|[2006l ) 
and the reddening law of lCardelli. Clavton fc MathisI (|l989t ) , 
we have obtained F(Ha) = 1.40 x 10~^^ ergcm"^ s~^ within 
0'.'4 of the nucleus. At the distance of 170 Mpc, L(Ha) = 
4.84 X lO^^ergs"^ and the corresponding nuclear bolomet- 
ric luminosity is Lboi = 4.84 x 10*'' ergs~^. If we assume an 
efficiency rj ~ 0.1, f or an optically thick and ge ometrically 
thin accretion disc (|Frank. King fc Raind 1200 2| ) we derive 
an accretion rate of rh = 8.7 x 10~^ yr~^. Since Mout is 
2 orders of magnitude higher than rh, we conclude that the 
outfiowing gas does not originate only from the AGN, but 
is composed mostly of interstellar gas from the surround- 
ing region of the galaxy, which is swept away by the AGN 
outfiow. 

We can additionally estimate the kinetic power of the 
outfiow, considering both the radial and turbulent compo- 
nent of the the gas motion. Using the H/3 kinematics to 
estimate this power, we have 



£0 



Mout. 2 
— ;^(*^out 



+ 0)., 



(13) 



iQyr % 
-5.0x lO'^'ergs'S 



where Dout is the velocity of the outfiowing gas and a is the 
velocity dispersion. From Fig.[5]we have a ~ 170 kms^^ and 
adopting t^out = 340kms~^ and Aiout = 0.3 — 1.1 Mr; 
as discussed above, we obtain iJout = 1.4- 
which is 0.03% - 0.1% x Lboi, so that between 0.03% 
and 0.1% of the mass accretion rate is transformed in ki- 
netic power in the outflow. These values are one order 
of magnitude smaller than the AGN feedback derived by 
iDi Matteo et al.l (|2005h in simulations to account for the co- 
evolution of black holes and galaxies and approxima t ely at 
the same order of those estimated by lBarbosa et al.l (|2009l ) 
and lStorchi-Bergmann et al.l (|2010l ) using similar IFU data. 
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6.6 Stellar Population 

The stellar population synthesis results are affected by the 
restrict wavelength range of our data, which excludes the 
near-UV region (below 4200 A) - the most sensitive to the 
intermediate age population. This limitation may have af- 
fected the absolute value of the contribution of this age com- 
ponent. In order to compensate for this, we have restricted 
the population vector to only 3 stellar components (young, 
intermediate age and old). The relative contribution of these 
different age components at the different locations seems to 
be robust, supported, for example, by the dominance of the 
young stellar population component in region 7, for which 
the emission-line spectrum supports ionization by stars. 

From the synthesis, we conclude that within the in- 
ner 0.3 kpc the old population dominates the flux at 4700 A 
(~60%), but that there is also some contribution of the in- 
termediate age one (~30%). The gaseous kinematics has 
shown that the AGN feedback is concentrated approxi- 
mately within the inner 0.3 kpc, coinciding with the region 
where the stellar population synthesis has revealed predom- 
inance of old population, but with some contribution of in- 
termediate age population and FC. Beyond this region, at 
the 0.8 kpc ring, the young stellar population has a contri- 
bution ranging from ~40% to «60%, reaching 100% at the 
detached region 7, located in an spiral arm. In regions 4, 5 
and 6, there is also significant contribution of the interme- 
diate age population. 

We can represent the mixture of ages of the stellar pop- 
ulation components in a more condensed form b y its mean 
stellar age. FoUowing lCid Fernandes et al.l (|2005l ) we can de- 
fine the mean stellar age weighted by the light fraction as 

N, 

< logt. >L= ^a^jlogij (14) 

j = i 

and by the stellar mass fraction as 

iV, 

< log >M= ^ mj log . (15) 
i=i 

While the light-weighted mean stellar age is biased to- 
wards younger ages, as we can see in Table [5] the mass- 
weighted stellar age is more representative of the older pop- 
ulation (wluch_Jws_a^Jajger ims^ ratio). Accord- 
ing to ICid Fernandes et al.l (|2005l ). although the latter is 
more physically representative, it has less relation with the 
observed spectrum. These authors argue that, in practice, 
< logf* >L is the most useful of the two indices, due to the 
large M/L ratio of old stars. In Table [2] we list these two 
mean ages for the integrated spectrum and for the regions 
1 — 7. < logt, >L is indeed characteristic of a post-starburst 
population at most locations, where mean ages are smaller 
than 1 Gyr, except at the nucleus, where the mean age is 
larger than 3 Gyr. 



7 SUMMARY AND CONCLUSIONS 

We have used integral field optical spectroscopy in order 
to map the stellar population and the gas kinematics within 
the inner 1.5 kpc (radius) around the Post-Starburst Quasar 
J0212-0903 at spatial resolution of ~ 0.5 kpc (0.6 arcsec) and 



velocity resolution of ~ 120 km s~ . This is the first two- 
dimensional study of both the stellar population and gas 
kinematics of a Post-Starburst Quasar and our main con- 
clusions are: 

(i) The stellar population is dominated by old stars within 
0.3 kpc (radius) from the nucleus, while in the circumnuclear 
region both intermediate age (lOOMyrs ^ t ^ 2.5G Gyrs) 
and young stars {t < 100 Myrs) dominate the optical flux. 
Our results support a location of most of the post-starburst 
population in a ring with radius of ~ 0.8 kpc, where active 
star formation is also occurring; 

(ii) Extended emission up to 1.5 kpc from the nucleus is 
observed in both [O iii] and H/3. The [O iii] flux distribution 
is more centrally concentrated, while the H/? flux distribu- 
tion is more extended approximately along the line of the 
nodes and is also observed in a detached patch which seems 
to belong to a spiral arm with active star formation; 

(iii) The emission-line ratios are typical of Seyfert excita- 
tion within the inner 0.3 kpc, where an outflow is observed, 
while beyond this region the line ratios are typical of LIN- 
ERs, which we attribute to a combination of diluted nuclear 
radiation and ionization by young stars. A Principal Com- 
ponent Analysis (PGA) supports the conclusion that recent 
star formation dominates the spectrum in the circumnuclear 
region; 

(iv) The mass of ionized gas in the inner 1.25 kpc (radius) 
is 8.34 ±0.92 x 10^ M©; 

(v) The gas kinematics can be reproduced by a combi- 
nation of (1) rotation in the plane of the galaxy whose line 
of the nodes runs approximately North-South with ampli- 
tude ^ 200 km s^^ (supported also by the PGA) and (2) an 
outflow observed within the inner 0.3 kpc, with blueshifts 
reaching up to — 670kms~^; 

(vi) The mass outflow rate is in the range of 0.3- 
l.lMeyr"^, which is ~ 100 times the AGN mass accre- 
tion rate of ~ 8.7 x 10"'^ M© yr~^, implying that most of 
the outflow originates via mass-loading in the surrounding 
interstellar medium of the galaxy, swept away by the AGN 
outflow. 

We conclude that our observations support both the 
evolutionary and quenching scenario for this galaxy, as fol- 
lows. The feeding of gas to the nuclear region has triggered 
a circumnuclear starburst a few lOO's Myr ago, extending all 
the way to the nucleus. The remaining gas from this inflow, 
combined with the mass loss from the newly formed stars 
may have then triggered the nuclear activity, producing the 
observed gas outflow. This outflow, observed within the in- 
ner 0.3 kpc, has then quenched further star formation at this 
location, in agreement with the observed absence of young 
stars and the contribution of intermediate age stars in the 
inner 0.3 kpc. Beyond this region affected by the outflow, 
star-formation seems to be still active in a 0.8 kpc ring. The 
presence of a delay between the triggering of star formation 
and the nuclear activity is supported by the presence of the 
intermediate age stellar population at the nucleus. 
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